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ABSTRACT 

Gamma-Ray Burst (GRB) afterglows are well described by synchrotron emission from relativistic blast waves 
expanding into an external medium. The blast wave is believed to amplify the magnetic field and accelerate the 
electrons into a power law distribution of energies promptly behind the shock. These electrons then cool both 
adiabatically and by emitting synchrotron and inverse Compton radiation. The resulting spectra is known to 
consist several power law segments, which smoothly join at certain break frequencies. Here, we give a complete 
description of all possible spectra under those assumptions, and find that there are 5 possible regimes, depending 
on the ordering of the break frequencies. The flux density is calculated by integrating over the contributions 
from all the shocked region, using the Blandford McKee solution. This allows us to calculate more accurate 
expressions for the value of these break frequencies, and describe the shape of the spectral breaks around them. 
This also provides the shape of breaks in the light curves caused by the passage of a break frequency through the 
observed band. These new, more exact, estimates are different from more simple calculations by up to a factor of 
~ 70, and describe some new regimes which where previously ignored. 

Subject headings: radiation mechanisms: nonthermal — gamma rays: bursts — gamma rays: theory — Shockwaves 



1. INTRODUCTION 

In recent years, several dozens of gamma-ray burst (GRB) af- 
terglows have been observed, and data is accumulating rapidly. 
The quality of these observations is constantly improving. The 
study of afterglow emission has helped shed light on many im- 
portant aspects of the GRB phenomenon. The spectrum dur- 
ing the afterglow phase is well described by synchrotron emis- 
sion from a relativistic blast wave, and consists of several power 
law segments (PLSs), that join at several break frequencies (e.g. 
Sari, Piran & Narayan 1998). These break frequencies are the 
self absorption frequency, v a , below which the optical depth to 
synchrotron self absorption is larger than unity, i/ m , the typical 
synchrotron frequency of the minimal electron in the power- 
law, and v c , the synchrotron frequency of an electron whose 
cooling time equals the dynamical time of the system. Granot, 
Piran & Sari (2000) have then found that if v c < v m , the self ab- 
sorption frequency actually splits into two: v ac and v sa , where 
an optical depth of unity is produced by non-cooled electrons 
and all electrons respectively. Different possible orderings of 
these break frequencies result in five possible spectral regimes, 
as shown in figure |IJ 

The physical parameters of a burst may be deduced from fit- 
ting the observed broad band spectrum to the theoretical spec- 
trum. This has been done by Wijers & Galama (1999) for GRB 
970508 by fitting a broken power law theoretical spectrum. A 
detailed description of the shape of the spectrum allowed a 
more accurate determination of the self absorption frequency 
v sa and the peak frequency v m (Granot, Piran & Sari 1999b, 
hereafter GPS99b). A more accurate theoretical calculation of 
the break frequencies leads to a more accurate conversion from 
the observed spectrum to the burst parameters. The combined 
effect was that the inferred value of the density, for example, 
was different than that of Wijers and Galama by two orders of 
magnitude. This illustrates the sensitivity of this method to the 
shape of the theoretical spectrum around the break points, and 



stresses the need for a more accurate determination of the theo- 
retical break frequencies, for all various spectral breaks. 

So far only the shape of the spectrum around v m (Granot, 
Piran & Sari 1999a, hereafter GPS99a; Gruzinov & Waxman 
1999) and v a (GPS99b) was calculate in detail, and even that 
was only done for the canonical case where v sa < v m < v c (see 
the upper panel of figure |l]). This paper, is intended to ex- 
tend these works for all spectral breaks, and therefore provide 
a comprehensive, self consistent calculation of the broad band 
spectrum. We provide analytic formulas which approximate 
the shape of each of the spectral breaks and their positions, in a 
form which is easy to use for afterglow fitting. We also suggest 
a prescription for combining these breaks to a single analytic 
broad band spectrum. 

The physical model is outlined in §|2| while a more detailed 
and formal description of the model and of the calculation of 
the observed flux density is given in §^ (Appendix A). Our 
main results are presented in §|^. In we give prescriptions for 
combining the shapes of the spectrum near the different spectral 
breaks into a single analytic broad band spectrum. We discuss 
our results in §g. 

2. THE PHYSICAL MODEL 

An exact calculation of the spectrum requires the knowledge 
of (i) The hydrodynamic quantities (bulk Lorentz factor and 
number density) (ii) The magnetic field strength (iii) The elec- 
tron energy distribution. These should be given for any location 
behind the shock, and at any time. Below we describe our ap- 
proach to all three. 

(i) The hydrodynamics is described by the Blandford Mc- 
Kee (1976, BM hereafter) self similar solution. This solution 
describes a spherical relativistic blast wave expanding into a 
cold medium, and assumes an adiabatic flow, i.e. that radiation 
losses are small and do not effect the hydrodynamics. Radia- 
tive effects can be taken into account to modify the hydrody- 
namic evolution as described by (Sari 1997, Cohen, Piran & 
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Sari 1998) and to modify the structure of the cooling layer be- 
hind the shock as described by Granot & Konigl (2001). If the 
radiative losses are not too large, our formalism would give the 
correct break frequencies and break shapes, provided that one 
uses the time dependent energy, as discussed in the first two of 
these references. The BM solution we use is for an impulsive 
explosion in an ambient density described by a power law with 
radius, p ext (r) = Ar~ k . We consider two different values of k 
which are of particular physical interest: k = 0, corresponding 
to an interstellar medium (ISM), and k = 2, corresponding to a 
massive star progenitor, surrounded by its pre-explosion wind. 
The assumption of a spherical flow is also adequate for a jetted 
flow at sufficiently early times, when the Lorentz factor of the 
flow is still larger than the inverse opening angle of the jet. We 
therefore have the complete hydrodynamic description in term 
of the total energy E, and the external number density n ext (or 
A in the case of wind). The hydrodynamic profile that is used, 
is given in equations [l4] through |l6[ 

(ii) We assume the magnetic field gets a fixed fraction, es, 
of the internal energy everywhere behind the shock, as given 
by equation [l2|. This would be the case if the shock amplified, 
randomly oriented, magnetic field decreases due to adiabatic 
expansion. Different assumptions on the evolution and orien- 
tation of the magnetic field were shown to have only a small 
effect on the resulting spectrum (GPS99a,b). 

(iii) The electrons are assumed to acquire a power law dis- 
tribution of energies, N(j) oc y~ p for 7 > 7 m i n ,o, immediately 
behind the shock. Their total energy immediately behind the 
shock is a fraction e e of the internal energy. After being accel- 
erated by the shock, the electrons cool due to radiative losses 
and adiabatic cooling. The former can be calculated from syn- 
chrotron theory and the latter from the density profile given in 
by the BM solution. The resulting distribution is given in equa- 
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Given the above, the observed flux density may be calculated 
as described in §^ (Appendix A). The spectrum for the optically 
thin breaks may be calculated using equation^, while equation 
34 applies more generally. 



3. RESULTS 

power law segments. All five possible different spectra, which 
are shown in Figure |l[ consist of between three and five differ- 
ent power law segments (PLSs). At each asymptotic PLS (suf- 
ficiently far from the break frequencies) we have F v oc v^t a . 
Altogether there are eight different PLSs, labeled A through H, 
from high to low values of (3. Note that there are two differ- 
ent PLSs with a slope of (3 = 1 /3. Both are produced by the 
low energy tail of synchrotron radiation, but in region D it is 
the non-cooled electrons that are responsible for the radiation 
while in region E it is the cooled electrons. Most PLSs appear 
in more than one of the five possible spectra (see Figure [j]). If 
one is only interested in the spectrum far enough from the break 
frequencies, then the normalizations of the different PLSs is all 
that is needed to accurately describe the spectrum. This is given 
in Table |[ The coefficients for PLSs A, G and H, slightly de- 
pend on p in a non analytic way. These were calculated for 
p = 2.2,2.5, 3, and a linear function (or a linear function multi- 
plied by an exponent) was used to describe the results. 

For PLS E, we find that the emission becomes dominated 
by the contribution from small radii (i.e. early times when the 
radius of the shock was small) for £>23/13~l .769. The elec- 
trons responsible for the emission in this regime have suffered 
considerable adiabatic cooling (as well as radiative losses). In 



this regime (k > 23/13) PLS E splits into two different PLSs, 
whose spectral slope (3 depends on k. Furthermore, the effec- 
tive size of the afterglow image, at a given observed time, in 
this regime depends on the observed frequency. However, since 
this new regime is somewhat out of the main stream of this 
paper, and in order to avoid confusion, we leave the detailed 
description of this new regime to a future work (Granot & Sari, 
in preparation). The normalization of PLS E, and the expres- 
sions for the spectral breaks b = 10, 11 (that involve PLS E), are 
therefor left out of Tables [j] and 0, respectively, for k = 2. 
break frequencies. The different possible combinations of the 
eight PLSs, result in 11 different break frequencies, labeled 
b = 1, 11 (see Figure |l]). Again, the same break frequency 
may appear in more than one spectrum. The values of the break 
frequencies, Vb, and the corresponding extrapolated flux densi- 
ties, F Vb ,ext, are defined at the point where the asymptotic PLSs 
meet. These can be directly calculated from the normalization 
of the PLSs that are given in Table [j], but for completeness 
they are given explicitly in Table |J The fit for the p depen- 
dence was redone in this table (with either a linear fit, an ex- 
ponent or a combination of the two) using the accurate results 
for p = 2.2,2.5,3, therefore resulting in slight (a few percent) 
inconsistencies with the previous table. 

shape of breaks. The Flux density near a spectral break, Vb, 
may be approximated by 
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where (3\ and are the asymptotic spectral slopes below and 
above the break, respectively, and s is a parameter which de- 
scribes the sharpness of each break. The sign of s is equal to 
that of (3\-f32 (i-e. positive if the spectral slope decreases across 
the break), while \s\ represents the sharpness of the break (the 
sharper the break, the larger |s|). The shape of most spectral 
breaks (except for b = 1,8, 10, 11), depends on the value of p, 
and so does the corresponding value of s = s(p), which is given 
in Table 2. All quantities which depend on the value of p were 
calculated for p = 2.2,2.5, 3, and are given in a form which is as 
exact as the functional parameterization permits at these values 
of p, and interpolates (or extrapolates) for other values of p, 
and should therefore be reasonably accurate for 1.5 <^p <^3.5 
(see discussion below equation [Tl| for p < 2). 

The break b = 1 has been investigated in detail by GPS99a, 
for k = 0, and they found that the physically motivated formula 



F v = F Vl ,ext [ 1 - exp(-Ti )] r, 6/5 , n = (v/vi ) 5/3 



(2) 



provides an even better description of F v near the break (with 
an MPvD of 2.63%, compared to 6.78% with equation|l]). How- 
ever, for k = 2, equation [l] provides a better fit (with an MRD 
of 1.02% compared to 25% with equation |J), which shows that 
the previous success of equation || was accidental. 

Equations [j] and |2| both give a poor fit for b = 4. This is 
because the spectral slope across this break does not change 
monotonically. We therefore provide an alternative formula for 
this break, 



F v = /v ext 4>l exp(-s</>4 /3 ) + (f> 5 / 2 



(j) 4 = V/V 4 



(3) 



where the values of s for b = 4 which appear in Table 2, are for 
this equation, rather than for equation [if, as for the other breaks. 
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4. A PRESCRIPTION FOR THE BROAD BAND SPECTRA 



5. DISCUSSION 



The values and the shape of the break frequencies, as given 
in the previous section, are strictly valid only when the break 
frequencies are far away from each other. Though, in principle, 
our formalism is adequate to describe the general spectrum, for 
arbitrary values of the break frequencies, such a description, 
would require a new calculation for any ratio of the break fre- 
quencies, and is therefore not practical. Instead, we choose to 
give a heuristic prescription, that uses the shapes from the previ- 
ous section to construct a broad band spectrum, which includes 
all the breaks, for an arbitrary ratio of the break frequencies. 
Once again we stress that this is not a rigorous derivation of 
such a spectrum, but simply an analytic equation, that gives a 
smooth spectrum when the break frequencies are close to each 
other, and approaches the rigorous shape of each break, in the 
asymptotic situation where the break frequencies are far apart. 
Such an equation is useful for fitting afterglow data. 

One can readily construct such a formula for any one of the 
five different possible spectra shown in Figure [j]. Let us label 
these spectra 1 through 5, from top to bottom, and denote the 
corresponding flux densities by F® where i= 1, ...,5. We also 
label the flux density near the 1 1 spectral breaks by F,, where 
b = 1 , . . . , 1 1 . The fluxes, F/,, are given by equation [j] (for break 
b = 4 equation || gives a more accurate description). Now, let us 
define a quantity Ft,, by 



F b = l + (v/v b 



-i/„ 



(4) 



The formulas for the rounded shape of the spectrum for the five 
spectra which are shown in Figure [jj, from top to bottom, are 
given by 



F$> = F,F 2 F 3 , 


(5) 


= F 4 F 5 F 3 , 


(6) 


F^=F 4 F 6 , 


(7) 


Fl 4) = F 7 F 8 F 9 , 


(8) 


F v = F7F10F1F9 . 


(9) 



The first term, F&, provides the normalization and the shape of 
the spectrum near the lowest break frequency, while each con- 
secutive term, F,, represents the next break frequency, from low 
to high frequencies, and provides the shape of the spectrum near 
that break frequency, and the appropriate change in the spectral 
slope across the break. The number of free parameters in each 
spectrum generally equals the number of break frequencies plus 
two, since besides the values of the break frequencies, one has 
to specify the value of p and of the flux normalization. The bot- 
tom panel of figure [j] is an exception, and F^ 5) has only 5 free 
parameters, since there is a closure relation between the four 
break frequencies (Granot, Piran & Sari 2000): 



^10 / v\\ 

Vl V ^9 



4/5 



= Hp) ~ 1 (for k = 0) 



(10) 



We have used the BM solution, to obtain more accurate ex- 
pressions for the flux density. Under the assumptions that the 
initial electron distribution is a strict power-law with a low en- 
ergy cutoff, and that the magnetic field is amplified immediately 
behind the shock, we derived exact expressions for the values 
of the break frequencies, as well as the shape of the spectrum 
around each break. We have given a complete general descrip- 
tion of the broad band spectrum. As our analysis is general, it 
also includes exotic spectra, that may only be relevant in very 
early phases or for extreme parameters. Our main results are 
summarized in Figure |l| and Tables 1 and 2. 

In general, the spectrum of GRB afterglows evolves from fast 
to slow cooling 3 . For example, for an ISM with standard pa- 
rameters, (e.g. «o — LF52 = 1, es — 0.01) the initial spectrum 
is 5, then v m crosses v c and the spectrum turns into spectrum 1, 
and finally, when v m crosses v sa the spectrum turns into spec- 
trum 2. The transition times between the various spectra of Fig- 
ure [l] can be worked out by equating the various frequencies, as 
given in Table || It follows that there are two types of evolution 
depending on the parameters, as given in the chart below. 



Our prescription for constructing the broad band spectrum 
uses F Vb ext from only one of the break frequencies in each spec- 
trum, and thus avoids the problem of the slight inconsistencies 
within Table || (that arise due to the independent fits for the p 
dependence, e.g. §||). 

3 This holds for k < 4, which includes the cases relevant for the afterglow, k = 0,2. 



ISM = 



WIND: 




4/79/7 

52 t B 

4/7 9/7 
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'Ft n a 1 > 100 
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-,-3/7 2/7 



eJ' < 100 



The weakest link in our formalism, is the assumption of a 
sharp lower cutoff in the electron distribution. This would affect 
breaks b = 1,2,4,7,9 (though for b = 1, the shape of the break 
will not be effected, while v\ and F, i ext may change). Nev- 
ertheless, our calculation provides the first self consistent de- 
scription of all these breaks. The values and shape of the rest of 
the breaks, depends only on the assumption of a power-law dis- 
tribution, well above the low energy cutoff, and on the electron 
cooling. Our description of these breaks (b = 3,5,6,8, 10, 11) is 
therefore more robust. These breaks may still be somewhat af- 
fected by the assumption of the magnetic field evolution. How- 
ever, in previous papers (GPS99a,b), we have shown that this 
dependence is relatively week (F, 2 . ex t and change by up to 
~ 50%, while F Vi >e xt and v\ change only by up to a few percent, 
where in both cases the shape of the break does not change con- 
siderably). 

We give a complete description of all possible power-law 
segments (PLSs), and provide exact expressions for the flux 
density away from the break frequencies. These expression are 
useful when partial information for the afterglow exists. In gen- 
eral a spectral slope and a flux at some frequency are sufficient 
to give some constraint on the afterglow parameters (in PLSs 
G and H, p would also be needed). For example, if only X-ray 
data exists, PLS H can be used to extract some information on 
the underlying parameters, and if only radio data exists, PLS 
B can be used, even if the self absorption frequency is not ob- 
served, (i.e. is above the observed radio frequency). 

Expressions for some of the break frequencies and corre- 
sponding flux densities already exist in the literature (Waxman 
1997; Sari, Piran & Narayan 1998; Wijers & Galama 1999; 
GPS99a,b; Granot, Piran & Sari 2000; Chevalier & Li 2000; 
Panaitescu & Kumar 2000). Most of these works address the 
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spectrum shown in the upper panel of Figure 1 (spectrum 1). 
The values we obtain for the break frequencies and correspond- 
ing flux densities are in some cases significantly different than 
previous estimates (by up to a factor of ~ 70). For k = and 
p = 2.5, our value for V2 (F U2 , ext ), which is better known as v m 
(F Vm ), is a factor of 1.3 (4.2) larger (smaller) than Sari, Piran 
& Narayan (1998), a factor of 1.5 (3.1) smaller (larger) than 
Wijers & Galama (1999), a factor of 3.8 (for vj) smaller than 
Panaitescu & Kumar, and a factor of 15 (8) smaller (larger) than 
Waxman (1997). For p = 2.2 our value for V2 is a factor of 70 
smaller than Waxman (1997) 4 . Our values for vi and F^ ext are 
only slightly different (by -5.1% and +1.6%, respectively) than 
GPS99a, due to a small approximation they made for the local 
emissivity. Our value for vt, {v c ) is a factor of 2.6 larger than 
Sari, Piran & Narayan (1998), a factor of 6.4 larger than Wijers 
& Galama (1999), and a factor of 6.1 larger than Panaitescu 
& Kumar. Our value for v\ {v sa ) is a factor of 1.9 larger than 
Waxman (1997), a factor of 3.7 smaller than Wijers & Galama 
(1999), a factor of 4.1 larger than Panaitescu & Kumar, and a 
factor of 2. 1 smaller than GPS99b 5 . For k = 2 (and p = 2.5), our 
values of v\, and ^3 are smaller by factors of 2.5, 1.4 and 
1.4, respectively, compared to Chevalier & Li (2000), while 
our value for F„ 2 ext is larger by a factor of 3.2; Compared to 
Panaitescu & Kumar, or value for v\ (y-i) is larger (smaller) by a 
factor of 4.9 (7.4), while our value for z/3 is larger by only 16%. 
Our expressions for the break frequencies and corresponding 
flux densities of spectra 4 and 5 (bottom two panels of Figure 
1), for k = 0,2, are different by up to a factor of 3 from those 
given in Granot, Piran & Sari (2000). 

Our equations do not include the effects of inverse compton 
scattering on the cooling of the electrons. This effect is known 
to be important, when eg <C e e (Sari, Narayan & Piran 1996, 
Panaitescu & Kumar 2000, Sari & Esin 2001). Following the 
prescription of Sari & Esin, we can include the effects of inverse 
compton by inserting appropriate powers of (1 +Y) into the val- 
ues of the break frequencies, or the PLSs (where Y is the comp- 
ton y-parameter). PLSs C, E, F and H should be multiplied by 
(l+Y)~y & , (1 +F) 2 / 3 , (1 +YT 1 and (1 +Y)'\ respectively. 

Preliminary results from this work have already been used 
successfully in fitting the data of several afterglows (e.g. 
Galama et al. 2000; Harrison et al. 2001). In the latter case 
the first evidence for inverse Compton emission was found. A 
special effort has been made to present the results of our model 
in a way that is simple to implement, and would provide the 
most accurate results to date for spherical afterglows, or jetted 
afterglows within their quasi spherical phase (before any signif- 
icant lateral spreading). 

JG thanks the Horowitz foundation for support. RS thanks 
the Sherman Fairchild foundation for support. This research 
was partially supported by a NASA ATP grant. 

6. APPENDIX A 



frame, i.e. the rest frame of the ambient medium in which the 
flow is spherical. We use a spherical coordinate system in this 
rest frame, where the z-axis points at the observer. The time, t , 
measured in this rest frame is called the coordinate time, and is 
to be distinguished from the time, t ', measured in the local rest 
frame of the fluid, and from the observer (or observed) time, 
f b S , at which the emitted photons reach the observer. The sub- 
script '0' denotes the value of a quantity just behind the shock. 

The initial electron distribution, just behind the shock, is 
given by 



N(-j e )=K j/ for 7, >7, 



£ee ° Ml, 
min.O = X , (H) 

riQm e c z 



where m e is the electron rest mass and Kq = (p— 1 ^olio- 
Note, that the above equation is usually written using e e = 
e e {p-\)/(p-2), which is the fraction of the internal energy 
given to the electrons. The advantage of using I e , is that it 
makes most equations somewhat simpler. Furthermore, it will 
apply also for the case p < 2, as long as the minimal Lorentz 
factor is proportional to the shock Lorentz factor. The magnetic 
field is assumed to hold a constant fraction, eg, of the internal 
energy, everywhere, 

B 2 = 87re B e. (12) 

The evolution of the Lorentz factor of each electron is de- 
scribed by 

dj e <J T B 2 2 , le dn 

~dl' 



6irm e c^ e 3n dt' 



(13) 



The first term on the right hand side of equation [13] represents 
the radiative losses while the second term represents adiabatic 
cooling. The radiative term includes only synchrotron losses. 
A simple prescription of how to include the effects of enhanced 
electron cooling, due to inverse compton scattering, on the ob- 
served synchrotron emission, is given in §B1. 

We use the BM spherical self-similar solution for an impul- 
sive explosion, where the external medium is cold and its den- 
sity changes as a power law of the distance from the center, 
Pext(f) =Ar~ k , k < 4 (extensions for k > 4 are given in Best and 
Sari 2000, but were not used in this paper). The derivations are 
made for a general value of k < 4, and are then used for k = 
and k = 2, which are of special physical interest. According 
to this solution, the proper energy density, Lorentz factor and 
proper number density of the shocked fluid are given by 



e = 2T p exl c x 

7 



-(17-4jt)/3(4-jt) 



= 2 -l/2p x -l/2 



« = 2 3 / 2 r« extX -< 1 °- 3 *>/ 2 < 4 - i >, 

where T is the Lorentz factor of the shock, and 



x =[i + 2(4-£)r 2 ](i-- 



The x coordinate of a fluid element is given by 



(14) 
(15) 
(16) 



(17) 



X 



4-k 



4-k 



(18) 



The energy density e, number density n, magnetic field B, 
and random Lorentz factor of the electrons j e , are measured in 
the local rest frame of the fluid, in addition to all the primed 
quantities. The remaining quantities are measured in the lab 

4 This large difference is mainly due to the fact that Waxman used e f instead of our e e = e e (p—2)/(p- 1). 

5 The reason for this last difference is as follows. Equation 18 of GPS99b, which is essentially equation 6.52 of Rybicki & Lightman (1979), misses the term 
associated with the discontinuity at the lower edge of the electron distribution (at 7 m ; n ) when derived from equation 6.50 of Rybicki & Lightman. This missing term 
caused an overestimation of the of the absorption coefficient by a factor off = 3(p+2)/4, and a corresponding overestimation of u\ and /Vj.exi, by factors of / 3 / 5 and 
/'/ 5 , respectively. However, this missing term does not effect the shape of the break, which is given in equations [l] or ^| (i.e. equation 24 of GPS99b). 
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where Ro and to are the shock radius and coordinate time, re- 
spectively, when the fluid element crosses the shock. Since 
T 2 oc t k ~ 3 , we obtain that 



7 (7-2A-) 

— — = y 2(i-t) 
70 



n <i3-2t) 

= V 2(4-4) 



e / B \ 2(13-2t) 

eo \BoJ X 

(19) 

Using e quat ion [T^ and the relation c/f ' = dt / 7 , we can write 
equationO in terms of \\ 



dx 



a T Blto X - (49 - m/6(4 - k) le 



6(4 — k)irm e cjo 
Solving equation [[o| we obtain 

"fe(je,Q,X) = 



(13-2*) 7e 
6(4-*) x 



7e,0 



x (13-2«/6(4-« + >i0/7rnax(x) 



(20) 



(21) 



where y e Q = %(x = 1) is the initial Lorentz factor of the elec- 
tron, just behind the shock, and 7 max (x) is the maximal Lorentz 
factor at x > 1, which corresponds to an electron with j e fi — ► 
oo, and is given by 



7max(x) : 



2(19 — 2*)7rm c c7o 



X 



(25-2*0/6(4-*:) 



^(19-2*0/3(4-*-) _ l 



(22) 



The fraction of electrons with a Lorentz factor within the inter- 
val [■j e ,j e + d'y e ] is given by: N( r y e )d / y e /n, and remains constant 
as all these quantities evolve with increasing %. The electron 
distribution is therefore given by: 



K qX 6<4-« % p l 



7. 



N(7e,X) 
p-2 



(23) 

for 7min (x) <le< 7max(x) 



7max(x) 

where 7 min (x) = 7e(7min,o, x)- 

We now have explicit expressions for both the hydrodynam- 
ical quantities and the electron distribution, over all relevant 
space-time, and can calculate the flux density near the vari- 
ous break frequencies. For breaks that are in the optically thin 
regime (b=2, 3,9,1 1) one may use the equation 



F v (tobs) = 



2{A-k)R]{\+ z ) 

d 2 



dy 
o Ji 



dx 



xy 2(5 - k) P'Ay,x,u s ) 

[l + (7-2fc)xy 



A-k] 



(24) 

which is a generalization of equation 13 of GPS99a, where 
o?£ and z are the luminosity distance and cosmological redshift 
of the source, respectively, P' v , is the radiated power per unit 
volume per unit frequency in the local rest frame of the fluid, 
and should be taken at the coordinate time t = t z + rfi/c, where 
f, = f obs /(l+z), 



H = 



(n-4k)(4-k)Et z 
AttAc 

(17-4*)£ 



i V(4-i) 



(25) 



1/2(4-*) 



_4 5 -*(4-£) 3 -*7rAc 5 -*f?-* 
E is the energy of the blast wave, y = R/R/ (e.g. GPS99a), 

1 — YV 4 ~ k 

fi = cos(6>) ^ 1 - xy 



(26) 



and 6 v' = 1/7(1— Pfj). The Spectral emissivity of a single elec- 
tron (in the fluid rest frame) is given by 



m e c 2 \u' sya 



3q e B r y 2 sin a 



(27) 

where q e is the electric charge of the electron, a is the pitch 
angle between the direction of the electron's velocity and the 
magnetic field, in the local rest frame of the fluid, and F is 
the standard synchrotron function (e.g. Rybicki & Lightman 
1979). In order to obtain an expression for P' v , (which appears 
in equation we average P' v , e over a, assuming an isotropic 
distribution of electrons in the local rest frame, 



,772 

K'eho = / damiaP' v , e {sma) 
Jo 



and then integrate over the electron distribution, 



d le N{-f e )Pl,^ (j e ) 



(28) 



(29) 



For the remaining spectral breaks (b=l,4,5,6,7,8,10), where 
the system is not always optically thin, we follow the formal- 
ism of GPS99b. Since the emission is isotropic in the lo- 
cal rest frame of the fluid, the emission coefficient is simply 
fv 1 = P'v' where P' v , is given by equation [29| The absorp- 
tion coefficient is given by 



1 



a„, = 



%-Ktn e v 



7™* N ( le) g 

a I rf T«— a7 be P v>,eMo(le)\ ■ (30) 



" mill 



Te die 



Since the flow is spherically symmetric, the afterglow image is 
circular, with physical radius of 



Ri 



_ (5-fc)zF*5 R, 
' ~ \/2 7/ " 

2 2 ~ k ( 1 7 - 4ifc)(4 - k) 5 - k Ec i - k t 5 - k 
TT{5-k) s - k A 



(31) 



1/2(4-*) 



and for a given observer time, t \,s, the specific intensity (or 
brightness), /„, depends only on the normalized radius from the 
center of the image, 



R 



R±,n 



■. (4 - ky l /2 (5 - *) ^ Vy-xy 5 ~ k 



(32) 



where x = at the center of the image and x = 1 at the outer edge 
of the image. As discussed in GPS99b, I v {x) may be obtained 
by solving the radiative transfer equation, 



dl u _ . _ T 
as 



(33) 



where s is the distance along the trajectory of a photon to the 
observer, and the flux density is given by 



4(4-*) 7/ 2 y ' 

6 Here, (3 is the fluid velocity in units of the speed of light, rather than the spectral 



F u (t obs ) = 2n(l+z) 

index. 



*R_L,max(?obs) 



dt 



xdxl„(x,t bs) ■ (34) 
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We note that I u (x) provides the surface brightness profile of 
the afterglow image, that is necessary for detailed calculations 
of microlensing or scintillation. The surface brightness profiles 
that were calculated according to this formalism, have already 
been used to study the microlensing of GRB afterglows (Granot 
& Loeb 2001; Gaudi, Granot & Loeb 2001), and are presented 



therein. 

When a break frequency is sufficiently far from other break 
frequencies, the spectrum near this break frequency assumes a 
self similar form. These self similar forms of the spectrum near 
the different break frequencies are presented in §|3. 
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Table 1 

The Normalization of the Different Power Law Segments 



PLS 


P 


F y (k — 0) in mjy 


F v (k = 2) in mJy 


A 


5/2 


l. 1 8(4.59- P )10«(l + ^ eB V4„-V2 £ i/4 f V4^ 8 ,5/ 2 


2.96(4.59-rtl0^1 + z)V4e-V^4/ 4 ^ s ^ 8 ^f 


B 


2 


4.20l|±fl 10»(1 + zfl^n ^ t^d^ul. 


1.33 10 9 (1 +z) 2 e ( ,A; 1 £52^^28^14 


C 


11/8 


8 oi ■ io 5 n +z) 27 / 1 V 1/ V 5/1 V /1 V 1/1 V 2 v n/ * 


3.28 ■ 10 5 (1 V' V'*^ V'&ViT 


D 


1/3 


27 9i£di(l + ,)V6 g -2/3 e l/3„l/2 £ 5/6 l/2 d -, 1/3 
5 \3p-l) w ^ z; te e fl "u c 52 'days a L28^14 


21li£dI(l +? )4/3 r 2/3 1/3^ £-1/3^-2 „V1 
Qp-iy- z > e B /1 * E 52 "L28^14 


E 


1/3 


73 On +z) 7/6 e n /6 E V6 t l/6 d~ 2 u l/i 
/J.U(.l-t-ZJ t B n Q -c 52 'days"L28^14 


... t 


F 


-1/2 


6 87H +^3/4 f -!/4 F 3/4 -1/4 ,_2 -1/2 
6.87(1 + z) eg t 52 t &iys d L2g u u 


6 fiSn 4-^V4 f -'/4 F 3/4 -1/4 ,_ 2 -1/2 
6.68(1 +z)'e B t 52 r days d L1% v u 


G 


(1-/0/2 


().461(p-0.04)f 2 - 53 "(l+z) — er e/ «J'\ 4 ? day 4 s d^u 


5+p l+fl l+p l-p 

3.82(p-0.18) e 2 - 54 "(l+z) — er'e B 4 A t E^ / days ^ 2 V, 4 2 


H 


-p/2 


0.855(p - 0.98)e L95 " (1 + z) T eV Ej 


2+1, I'-- 1 2 +l> 2 -'!' n 
0.0381(7.11-p)e 2 - 76 "(l+z) — ef V E 5 * df^v* 1 



NOTE. — The first two columns give the labels and the spectral slope, j3, of the different PLSs (see Figure MJ), while the last two columns give the asymptotic flux 
density within each PLS, for k = and k = 2. The reader is reminded that e e = e e (p—2)/(p— 1) depends on p. Trie notation Q x stands for the quantity Q in units of 10 r 
times the (c.g.s) units of Q, while ^ays i s the observed time in days, and A* is A in units of 5 • 10 gr/cm (Chevalier & Li 2000). 

t For PLS E,the emission becomes dominated by the contribution from small radii for k > 23/13. This new regime is described in a separate work (Granot & Sari, 
in preparation). 
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FIG. 1 . — The different possible broad band synchrotron spectra from a relativistic blast wave, that accelerates the electrons to a power law distribution of energies. 
The thin solid line shows the asymptotic power law segments (PLSs), and their points of intersection, where the break frequencies, i%, and the corresponding flux 
densities, iv 4 . e xt, are defined. The different PLSs are labeled A through H, while the different break frequencies are labeled 1 through 11. The temporal scalings of 
the PLSs and the break frequencies, for an ISM (k = 0) or stellar wind (k = 2) environment, are indicated by the arrows. The thick solid line shows the spectrum we 
calculated in this paper, where the broad band spectrum is constructed according to the prescription suggested in 84. The different spectra are labeled 1 through 5, 
from top to bottom. The relevant spectrum is determined by the ordering of the break frequencies. The top two panels (spectra 1 and 2) correspond to slow cooling 
(u,„ < u c ). Spectrum 1 applies when v sa < v m , while spectrum 2 applies when v m < v sa < v c . The two bottom panels (spectra 4 and 5) correspond to fast cooling 
(u c < u,„). Spectrum 5 applies when v sa < u c , and spectrum 4 applies when v c < u m < u,„. Spectrum 3 (middle panel) applies when v sa > v m ,v c , where in this 
case the relative ordering of v c and u,„ is unimportant (i.e. spectrum 3 may apply both to slow cooling or fast cooling). 
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Table 2 

The Break Frequencies and Corresponding Flux Densities 



b ft fa v b 



Vb(p) in Hz 



Fu b: ea(p) in niJy 



■<P) 



MRD in % 



1 21 (P ~' )V5 10 9 (1 I .-r'c-'c l/5 n 3/5 £' /5 
3 Vs " £31 ( "-" 3/5 lQ9fl | rr 2/5 c-'c' /5 .l 6/5 £" 2/5 r 3/5 

(3p+2) 3 / 5 1 U l 1+ *-f £ f e B /1 * C 52 Mays 



3.73(p-0.67)10i5(l + z)i/2 £ V^2 e V^ 

1/2-2 1/2,-3/2 



Q.647 (p -' )6/5 l/5 (l + Z ) 1 / 2 e -- 1 ^« /5 «n / ' 4 / '°'H 1/2 

Qp-llO/H-l) 1 / 5 k e s 52 days L2: 

9.19 (p -' )6/5 1/5 (l +Z) 6 '%- 1 el^Ell^dT* 
(3p-l)(3p+2)'/ 5 V e s * 52 days Z.28 



9.93(p+0.14)(l + Z )e7 1 ^ 72 £ 52 rf-f 8 
76.9( P+ 0.12)(l +z )V 2 ^ 2 A t £ 5 y fda 1 y / 2 rf Z | 8 



1.64 
1.06 



6.68 
1.02 



1 1-p 



4.02(p-0.69)10 l5 (l+ Z ) 1/2 £'/ ! e]e 



days 



,, +] j p-lE ^' j ~' J " 

n E S2 f days"L28 



1.84-0.40p 
1.76-0.38p 



5.9 
7.2 



1-P _P 



6.37(p-0.46)10 13 e-' 16 ' , (l+zr l / 2 £ B 3/2 n 1 £' 52 1/2 f da 1 y / s 2 
' c 4.40(3.45-p)10 10 e a45 "(l +zr 3 / 2 e B 3/2 A; 2 £ 5 ' 2 /2 4 /2 



4.68e 4 ' 82(,, - 2 - 5) 10 3 (l+z)T"e 
802g 7.02( P -2.5) 10 5 (1+z)P 4g-l £ P-i A ; £ l/2 ; i-P t/ -2 8 



e e tfl "0 c 52 'day 

3.04(p- 1.79)10 15 (1 +z) 3 ^e B A;'£ 2 f . 2 rf7 2 



1.15-().06p 
0.80-().03p 



1.9 
4.4 



4 2 



5 5WP- 1.22)10^(1 +Z) 1 / 2 .14 / 4 Xys 

2 " M 8.08(p-1.22)10 | 6(l+z)V2 e -24/2 £ l/2 f -3/2 



3.72(p-1.79)10 15 (l +z ) 7 / 2 6l6 fl n ''^\-;';j Z28 



3.44p-1.4lt 0.7t 
3.63p-1.60t 1.8+ 





' 4(„-l) p+2 4 £ p+2" 


l/2(p+4) 




'<i + c) 7 " +3 4" +3 4' +7 " 








3.59(4.03-p)10V 34 ' J 


e e s fl n E 52 




20.8(p-1.53)e 2 ' 56 'V7j 8 




1. 47-0.2 lp 


5.9 


(1+Z) 6 -Pt] p+2 
v ' days 




-10(l-p) 5(p-l) 
« days 





5 l-£ 
2 2 



1.58(4.10-p)10 lu e 



10-2,16p 



-i l/2(p+4) 



(i+j) 2 -"£;:''r 3(p+2) 

x ' 52 days 



158(p-1.48)e 2 ' 24 ''J7 2 



iu.i l l»', 1 '*V'" 1 

-10(l-p) A 2(;,-6) ( 4 ; ,+ l 
e * days 



1/20*4) 



1.25-0.18p 



7.2 



3.23(p-1.76)10 12 
4.51(p-1.73)10 12 



4(p-l) p-1 ,4 
days 



1725+5) 



l/2(p+5) 



76.9(p-1.08)e 2 06 ''rf7j 8 



78.6(p-1.12)e L89 P^f 8 



(i+c) 7 " +5 4"" 5 £ 5 3 r 5 

IW-,l,i(Hi 

« days 

<l+z) 6 '' +5 e 7 ''" 5 £ 4 ;' t5 
J0(l-p );J 2p f 4p-5 
days 



i/2(P+5) 



1/20*5) 



0.94-0. 14p 
1.04-0.16p 



e * days 

5 27 (3p -"" /i 10- 3 (1 I -rl/lo ? -" /5 c~ 4/5 » 1/10 £ 3/10 f" / "'i- 2 
(3p+2)»/5 £ f e B "0 C 52 'days "L28 

(3p+2)"/ 5 ' B A * C 52 f d-'ys"L28 



12.4 



11.0 



11 



1.6 8^gl O^(l + z)-'^ 5 ^ 5 A^ /5 



1.99-0.04p 
1.97-0.04p 



i.98.io"(i +z) - i / 2 »; /b £ 5 i rc 

8 2 "» 3.15-10"(l+zr l/3 Al /3 f da 2 y /3 



154(1 +z ) e B 1/ V /12 i^ 3 d7j 8 



1.9 
1.9 



Vm 3.52(p-0.31)10 1 5(l +z) l/2g2 e l/2 E l/2^ 



_____ 32 days 

0.221(6.27-p)(l+ Z ) 1 / 2 6- 1 efi ,/2 £ 5 1 2 /2 f da ,2 rf Z 2 8 
0.165(7.14-p)(l+z)V2 e --i e -i/2^2^^ 



0.907 
0.893 



1.71 

2.29 



3.34-0.82p 
3.68-0.89p 



7T~i i.32.io">(i +z r'/ 2 4 /5 »i 1/10 4 /1 \a 1 y / s 2 

8 3 v " ... tt 



3.72(1 + z )e 7/5 ^ /5 £ 5 7 2 /5 a f i ; 2 8 
... tt 



4.5 
4.2 
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1.213 
... tt 



5.86.io' 2 (i +z r'/ 2 , fi 3/2 V4 /2 C / " 

3 2 Vc ... tt 



28.4(1 + z )£ B /2 «y 2 £ 5 2^8 
... tt 



5.22 
... tt 



0.597 
... tt 



0.55 
... tt 



NOTE. — The first column numbers the breaks. The following two columns are the asymptotic spectral slopes below (/Ji) and above (/3 2 ) the break. The fourth 
column gives the name of the break frequency. The following two columns are u/,(p) and F Ub<: aip). The last two columns are the parameter s(p), which determines 
the shape of each break according to equation N (except for b = 4, where it applies to equation R|), and the maximal relative difference (MRD) between this analytic 
formula and our exact numerical results. For each break frequency there are two lines, the first is for an ISM surrounding (k = 0) and the second for a stellar wind 
environment (k = 2). The reader is reminded that l e = e e (p—2)/(p—V) depends on p. 

tFor b = 4, the values of s(p) and the corresponding MRD refer to equation and not to equation [l] as for the other breaks. 

ttThe breaks b = 10, 11 involve PLS E, where the emission is dominated by the contribution from small radii for k > 23/13. This new regime is described in a 
separate work (Granot & Sari, in preparation). 
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Table 3 

The Transition Times Between the Different Spectra 





possible definitions 


k 


transition time, ti-tj, in days 


5^ 1 


1/2 = vi, vg = vi\,vj = v\o 




2 


7.3- 10 2 -1.7- 10 3 x(l+z)e 2 e e 2 B n E S 2 
2.0 • 10 2 -7.0 ■ 10 2 X (1 +z)e e e B A* 


1 -> 2 


V\ = "2, ^4 = 




2 


6.1 • 10 4 - 1.2 ■ 10 6 X (l+ z )e 2 4 /5 n p 2/5 £' /5 

1.2 • 10 7 -3.9 ■ 10" X (1 +z)^ 0/3 ey W? 3 E 52 




i/ 10 = i/nt 


2 


9.3(1 +z)el /7 AlE^ n 


4^3 


f4 = ^6. ^7 = ^8 = ^9 




2 


2.2- 10 4 - 6.3 ■ 10 5 X (l+z)e 2 e e l B /2 n' 1/6 E 1 ' 3 
1.5 • 105-1.1 ■ 10 7 x (1 + z)el 2 /5 ef A^E% 5 


3^2 


^3 ="5 




2 


2p+l p+ti p+5 

5.1 ■ 10 8 - 1 .2 ■ 10" X (1 +z)e 2 e/ rr nf 1 

2(,,-l) 2p+7 2(p+6) 3 

8.0-24 x(l+z)e e 2 " +5 ef^A^ 5 E^ 5 



NOTE. — The first column indicates the transition at hand, from spectrum i to spectrum j. The second column lists possible conditions that may be used to define 
the transition time. The third column is k, which is either or 2, for an ISM or stellar wind environment, respectively. The last column is the transition time, ti-tj. 
There are several different ways to define most of most transition times (see second column), resulting in numerical coefficients that differ by a factor of order unity. 
The p dependence also varies the numerical coefficients by a factor of order unity. We specify the range of the numerical coefficients for 2.2 < p < 3 and for the 
different definitions of each transition. 

tThe expressions for v w and u n for k = 2, that we used in order to calculate t 4 ^ 5 , are taken from Granot & Sari (in preparation). 



